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Abstract 

We investigate the family of electrostatic spherically symmetric solutions of the 
five-dimensional Kaluza-Klein theory. Both charged and neutral cases are consid- 
ered. The analysis of the solutions, through their geometrical properties, reveals the 
existence of black holes, wormholes and naked singularities. A new class of regular 
solutions is identified. A monopole perturbation study of all these solutions is car- 
ried out, enabling us to prove analytically the stability of large classes of solutions. 
In particular, the black hole solutions are stable, while for the regular solutions the 
stability analysis leads to an eigenvalue problem. 
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1 Introduction 



Five-dimensional Kaluza-Klein theory [l], [J, or sourceless general relativity in 4+1 space- 
time dimensions (the extra space dimension being compactified) , was historically one of 
the first unified field theories. While this simplest higher-dimensional field certainly 
cannot be considered as realistic, it nevertheless deserves to be investigated as the pro- 
totype of other multidimensional theories. The static, spherically symmetric solutions 
of Kaluza-Klein theory have been obtained independently by several authors J3], f|, ||. 
These solutions include regular black holes, which are generalisations of Schwarzschild 
black holes with electric and scalar charges. A systematic investigation of these black 
hole solutions was carried out by Gibbons and Wiltshire 0. A class of regular, horizon- 
less charged solutions with wormhole spatial topology was also identified by Chodos and 
Detweiler in ||. Apart from these black hole and wormhole solutions, all other solutions 
apparently possess naked singularities. 

The aim of the present work is to analyze more fully the geometrical properties of the 
static spherically symmetric solutions of Kaluza-Klein theory, as well as to investigate 
their stability. As far as we know, the first systematic examination of the stability of 
these solutions is that of Tomimatsu J7|, which was restricted to electrically neutral solu- 
tions. Tomimatsu concluded that the only stable neutral solution was the 5-dimensional 
embedding of the Schwarzschild black hole. The stability of a class of wormhole solutions 
was analysed in ||; the conclusion was that these were generically stable. This investi- 
gation was generalized in an unpublished work |5| to encompass all the static spherically 
symmetric solutions. However, our subsequent analysis of the stability of scalar-tensor 
black holes |TIJ led us to uncover a flaw in the arguments of ||, || , which motivated us to 
launch a systematic reinvestigation of the Kaluza-Klein problem. 

In the second section of this paper, we recall the construction of the static, spherically 
symmetric (in the three "external" space dimensions) solutions to Kaluza-Klein theory. 
The solutions depend generally on three parameters x, a and b, related to the mass, 
scalar charge and electrical charge. The analysis of the geometrical invariants and the 
geodesies, carried out in section 3, allows the identification of the different solutions as 
black holes, wormholes and naked singularities in the five dimensional space-time. Some 
of these naked singularities turn out to be at an infinite geodesic distance, so that the 
corresponding solutions are regular. The equations for small time-dependent monopole 
perturbations of these solutions are then set up in a gauge-independent fashion in Sect. 
4, and decoupled in a special gauge. The transformation of these perturbations under 
gauge transformations is also briefly discussed. 

Using this framework, the question of the stability of static solutions to Kaluza- 
Klein theory is reduced to an eigenvalue problem. An analytical investigation of this 
problem is carried out in Sect. 5. We are able to prove that two classes of solutions are 
stable. The first stability class includes (contrary to Tomimatsu's claim) all the electrically 
neutral solutions, while the second stability class includes, among others, all black hole 
and extreme black hole solutions. We then discuss two special cases in which we are able 
to prove unstability (the second one containing a stable subcase). In the remaining cases, 
where we have no analytical information about the spectrum of eigenmodes, one should 
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resort to numerical computations to ascertain whether the corresponding solutions are 
stable under monopole perturbations. 



2 Dimensional reduction and electrostatic solutions 

The field equations of Kaluza-Klein theory derive from the 5-dimensional Einstein- 
Hilbert action 

S-- 



V [d 5 x^\R 5 , {2.1) 



16vr( 

with the additional assumption that d/dx 5 is a Killing vector with closed orbits. This 
last assumption allows the 5-dimensional metric to be decomposed as 

dsl = g„ u dx» dx v - e 2 ^{dx 5 + 2A fJ/ dx^f . (2.2) 

The 5-dimensional Einstein equations then reduce to the 4-dimensional system 

R„ v = -2e%F" + e-^^V, (2.3) 



D„ ( 



e^F^j = 0, (2.4) 
De V = -e^F^F^ , (2.5) 

with Fy, v = A Vi(t - A^ v . 

The first of these equations exhibits a non-minimal coupling of the scalar field e^ 
to 4-dimensional gravity. A minimal coupling is recovered by making the conformal 
transformation g^ u = e~^g^ u to the Einstein frame, leading to 

dsl = z~ 29l ^9% dx v - e 4 ^( dx 5 + 2A fl dx^f , (2.6) 

with the dilaton field tp = y/3tjj/2. This Einstein frame, frequently used in the literature 
0, is not defined when g 55 = — e 2 ^ is not negative definite, which will be the case for a 
large class of the solutions discussed in this paper. For this reason, we will use only the 
string frame defined in (|2.2|) . Let us also note that for A^ = the equations (|2.3|) 
reduce to the Jordan-frame equations of Brans-Dicke theory for uj = [11|. 



The spherically symmetric, electrostatic solutions of Kaluza-Klein theory have been 
previously obtained by several authors || [|, [5j. To be self-contained, we shall rederive 
them here along the lines of PJ, using the Maison approach to the dimensional reduc- 
tion of the higher-dimensional Einstein equations [|l^], which we first briefly summarize. 
The metric of (n + p)-dimensional spacetime with p commuting Killing vectors may be 
parametrized by 

ds 2 = X ab {dx a + ai dx l ){dx b + a) dx ] ) + r^hij dx l dx j , (2.7) 

where i — 1, n, a — n + 1, n + p, r = |det(A)|, and the various fields depend only 
on the coordinates x l . In our case n = 3, p = 2, x 4 is the time coordinate and x 5 is 
the Kaluza-Klein angular coordinate. Using the (n + p) -dimensional Einstein equations, 
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the magnetic-like vector potentials ft" may be dualized to the scalar twist potentials u a 
according to 

co a>i = |/i|-VVA^V ifc ^V ( 2 - 8 ) 
The remaining Einstein equations may then be written as the n-dimensional Einstein-cr- 
model system 

= 0, (2.9) 

Rij = \'R(x~ 1 X,iX~ 1 Xj) , (2-10) 

where the n-metric is hij, and x is the (anti-)unimodular symmetric (p + 1) x (p + 1) 
matrix-valued field 

*=(^CT'C"*)' < 2 - n > 

Solutions of Eq. ( [2.9|) depending on a single potential are geodesies 

X = Ve Aa (2-12) 

of the target space SL(p + l,R)/SO(p + 1), r\ and A being real constant matrices (with 
| det?7 1 = 1, Tr(v4) = 0, rj T = 7], A t t] = rjA) and a a harmonic function, 

VV = 0. (2.13) 

Now we specialize to n = 3, p = 2, and restrict ourselves to electrostatic solutions, 
uj a = 0. If cx(oo) = 0, the metric ( |2.7| ) is asymptotically Minkowskian provided 

/I 0\ 

77 = -1 . (2.14) 
\o \) 

Then the 3x3 matrix A is block-diagonal and may be parametrized by 

For a spherically symmetric solution with the potential a(r) normalized by lim( r ^ QO ) ra(r) = 
— 1, the parameters x, a and b are related to the physical observables M (ADM mass), £ 
(scalar charge) and Q (electric charge) by [0, [lj 



x = 2(M - S/v^) , a = -4S/v / 3, & = 2Q . (2.16) 
We may choose for the spherically symmetric reduced spatial metric the parametrization 

- dx i dx j = dr 2 + H(r) dtt 2 . (2.17) 
The harmonic function a(r) then solves 

o )T = H- x . (2.18) 
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Inserting the ansatz ( |2.17| ) in the reduced Einstein equation ( |2.10|) , we obtain the system, 

1 — -H >rr = , 

-H~ l H, rr + -H~ 2 H 2 = 2u 2 H~ 2 , (2.19) 
2 



with 
where 

This system is solved by 



v 2 = \ (x 2 - y) , (2.20) 

y = detN = b 2 + ax - a 2 . (2.21) 

H(r)=r 2 -v 2 . (2.22) 

The form of the function a(r) obtained by integrating Eq. ( [2.18 ) depends on the sign 
of the constant v 2 : 

a) y < x 2 (v 2 > 0). 
In this case, 

a= — ln(?— ^) (2.23) 
2v \r + v) y ' 

diverges for r = v. From the form of the function r = e xcr , we see that r = v is for x < 2u 
(x < or x > 0, y < 0) a point singularity of the 5-dimensional metric, and for x > 2u 
(x > 0, y > 0) a Killing horizon. 

b) y = x 2 (v 2 = 0). 

In this "extreme" case the function 

<r = -~ (2.24) 

r 

diverges for r = 0, which is a point singularity if x < 0. 

c) y > x 2 (v 2 < 0). 

In this case, which corresponds to the "Class III" of [El], the function 

^-Ig-^taQ), (2.25) 

where 

V 2 = = ^(iV-a: 2 ), 

varies within the finite range I a =] — ir/fi, [ when r varies from r — > — oo to r — > +oo. 
We therefore anticipate that the resulting 5-dimensional metric will be exempt from 
singularities. 

The full 5-dimensional metric (|2.7|) may be written down once the matrix A(r), given 
by (|2.12Q and ( |2.11|) , is known. The exponential in ( |2.12j ) is easily computed with the aid 
of the Lagrange interpolation formula |15| 

f(A)=J2f(z l )Y[^^- (2-26) 
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Its form depends on the sign of 

q 2 = x 2 /A-y, (2.27) 

which determines whether the eigenvalues of the matrix A are real, degenerate or 
complex. We obtain in these three cases: 



A 



1 



(fx 



X \ 

— — a J sinh qa + q cosh qa 
b sinh qa 



b sinh qa 
a ) sinh qa — q cosh qa 



for y < 



x 



4 ' 



A 



^Xcr/2 



X 



— a ) a + 1 



6a 



.r 



A = - e x,T / 2 
P 



— a ) sin pa + p cos per 
6 sin pa 



ba 



a I a — 1 



6 sin pa 



— — a j sm pa — p cos pa 



for |/ 



4 ' 



for -u > — , 
y 4 ' 

(2.28) 



where p 2 = y — x 2 /4. In all cases it can be checked that r = e 2 ^. We will also use, besides 
the 5-dimensional metric ( |2.7| ), the reduced 4-dimensional fields g^ u , and e 2 ^ obtained 
from (12.21), which we shall write in the form 



with 



/la 



ds 2 = e 27 dt 2 - e 2a dp 2 - e 2/3 dQ 2 , A^dof = V dt , (2.29) 



-e-^V* e 2/3 = #e-* CT , e 27 = e^ , e 2 ^ = -A 55 , V = A 45 /2A 55 . (2.30) 



If p in ( |2.29|) is identical with the radial coordinate r defined in ( 2.17]) , e 2a = e xa . 



3 Black holes and regular solutions 

3.1 Case y < x 2 . 

As observed in the preceding section, the surface r = v [a — ► — oo) where r = |det(A)| 
vanishes is a Killing horizon of the 5-dimensional metric, as well as of the reduced 4- 
dimensional metric ( [2 .29 ), for x > 0, y > 0. However the 5-dimensional metric ( |2.7| ) 



and the reduced 4-dimensional fields (|2.30|) are analytical only for y = 



On the Killing horizon r = z/, the matrix A reduces for y = 0, a; > to 

A* = 1 f " ) , (3.1) 



with eigenvalues and (2a — x)/x. If a < 0, then 2a < < x, so that the non- 
zero eigenvalue (2a — x)/x is negative and the non-null Killing directions are spacelike, 
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corresponding to an event horizon. The black-hole solutions of this class ||, which 
includes the 5-dimensional Schwarzschild metric for a = b = 0, have been extensively 
studied by Gibbons and Wiltshire ||. On the other hand, if a > 0, then (from y — 0) 
x = a — b 2 /a < a < 2a, so that the non-null Killing directions are timelike, and the 
Killing horizon is actually a conical singularity. The neutral solution 6 = (x = a) 

rfs 2 = dt 2 _ (cfe5)2 _ r + u_ dp 2 d[}2 

r + z/ v ' r -v v ' v ' 

is the direct product of the time axis by the Euclidean Schwarzschild metric, and is 
regular if the period of x 5 is 8ttu ||. The general b ^ solution behaves near r = v as 
the product of the (2+l)-dimensional metric generated by a spinning point particle |nj 
by the 2-sphere S2. 

In the case y 7^ 0, it is generally assumed that r = v is a curvature singularity [00. 
This assumption is based on an extrapolation from the neutral case 6 = 0. However 
the extrapolation goes through only if y < x 2 /4, in which case the matrix A has real 
eigenvalues and can always be diagonalized, so that the 5-dimensional curvature invariants 
of a charged solution are identical with those of some neutral solution. On the other hand, 
if y > x 2 /A, then the 5-dimensional metric cannot be diagonalized, and we must compute 
the metric invariants to ascertain whether they diverge on the Killing horizon. 

The components of the 5-dimensional curvature tensor may be computed from the 
formulas given in 0. Using the parametrization (|2.15| ), (|2.17 ) and Eq. ( |2.18|) , we obtain 

Rabcd = ~ tH- 2 [{XN) ac (XN) bd - (XN) ad (XN) bc ] , 

Ram, = \ H~ 2 [A (N 2 + (x - Ar)N)] ab , 

Raeeb =^H-\2r-x)(XN) ab , (3.3) 

Rr-ere = ~ T~ l H~ l {y - x 2 + 2xr) , 

R 6ip e v = 1 T" -1 sin 2 9 (2x 2 - y - Axr) . 

The computation of the 5-dimensional Kretschmann invariant K§ = R al3 " f6 R a ^s is easily 
carried out by tracing powers of the matrix N, for instance 

R ab cd R cd ab = I r 2 H-* ([Tr(iV 2 )] 2 - Tr(iV 4 )) , (3.4) 

leading to the compact result 

K 5 = 1 t 2 H~ 4 [y 2 + Apxy + 8p 2 (x 2 - y)\ , (3.5) 

with p = r — x/2. 

^From ( |3.5| ) we obtain the behaviour of the Kretschmann scalar near the Killing hori- 
zon (r ~ z/), 

_ 3(x - 2v) 2 (2x 2 - 3y) fr - v\ (*-**)/» 
5 ~ 512^ \~lh7 



(3.6) 
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For v > x/A (y < 3x 2 /A) this diverges on the Killing horizon, except in the special cases 
v = x/2, which corresponds to the y = black holes, or y = 2x 2 /3 f\ On the other 
hand, for v < x/A (3x 2 /4 < y < x 2 ) the 5-dimensional Kretschmann scalar is finite or 
vanishes on the Killing horizon. The form of the Kretschmann scalar K± for the reduced 
4-dimensional metric ( 2.29 ) is more involved as it depends on that of the Kaluza-Klein 



scalar field A55. We have found that on the Killing horizon r = u, it likewise diverges for 
y < 3x 2 /A (except for y — 0), and is finite or vanishes for 3x 2 /A < y < x 2 . 

Finite curvature invariants often, but not always, signal regularity of the spacetime 
geometry. A recent counter-example is that of "cold black hole" solutions [II], IT7), which 
have everywhere finite curvature invariants but, because of a lack of analyticity which 
prevents a Kruskal-like extension through the Killing horizon, are generically singular, 
except for a discrete set of solutions. In the present case we will see that the 3x 2 /A < y < 
x 2 solutions are indeed regular, although non-analytical on the horizon. They need not 
be extended through the Killing horizon because this is at infinite geodesic distance. In 
other words, they are already geodesically complete. 

The energy integral for geodesic motion in the spherically symmetric 5-dimensional 
metric ( p . 7|) - ( pTT7|) may be written as 



l2 2 

f2 + ^ + y W = ' v = - n T A- x n) , (3.7) 

where II4, II5 and I are constants of the motion proportional to the energy, electric charge 
and angular momentum of the test particle, r = dr/dp with p an affine parameter, and 
the (arbitrarily scaled) integration constant rj is positive, zero or negative for timelike, 
null or spacelike geodesies. Assuming x 2 /A < y < x 2 , x > (the Killing horizon can be 
shown to be at finite geodesic distance for y < x 2 /A, x > 0), we obtain from ( |2.28| ) the 
effective potential 

V = e xa (77 - ae- xa/2 sm(pa + /?)) (3.8) 

where the constants a and /3 depend on II4 and II5. 

For / = 0, e~ xa l 2 increases when when the horizon r = v (a — > —00) is approached so 
that, whatever the values of the constants of the motion, with a^O, there is a U\ such 
that 

e~ xai/2 sin(j9(Ti +/3) = rj/a , (3.9) 

corresponding to a turning point of the geodesic (reflection on a potential barrier). The 
only geodesies which can reach the horizon are those with a = (II4 = IK- = 0) and 
i] < 0, leading to 

p~ (r - i/)* 4 "-*)/ 4 " . (3.10) 

It then follows that: 

a) If x 2 /A < y < 3x 2 /A (x < Av), those geodesies terminate on the horizon (on which, 
as shown previously, the curvature invariant K 5 diverges). The 5-dimensional spacetime 
is geometrically singular, although this singularity is harmless to physical test particles 



5 However the fact that in this case the Kretschmann scalar vanishes on the horizon does not seem 
relevant, as the 5-dimensional metric is not analytical and cannot be extended through the horizon. 
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(V > 0), which will be reflected away before hitting the singularity (the same can be 
shown to be true for y = x 2 /A if x < 2a). 

b) If 3x 2 /4 < y < x 2 (x > 4i>), the horizon (on which Kf 5 \ is finite or vanishes) is at 
infinite geodesic distance, so that the 5-dimensional spacetime is geodesically complete. 
While these 5-dimensional spacetimes are geometrically regular, they admit closed time- 



like curves. From (|2.28|) , A55 changes sign periodically, so that all the circles r = 6 = 
ip = t = const, are timelike in the domains where A55 > 0. This is linked with the fact 
that the reduced 4-dimensional metric g^ u in (|2.2| ) is singular on the spheres A55 = 0, 
which separate 4-Minkowskian regions (A55 < 0) from 4-Euclidean regions (A55 > 0). We 
take the view that these 4-dimensional singularities are artefacts due to the breakdown 
of 5-to-4 dimensional reduction, i.e. the choice of a bad coordinate system (|2.2|) for the 
5-dimensional geometry, which is perfectly regular in the parametrization (|2.7|). 



3.2 Case y = x 2 . 

In this case v = and a = — 1/r. We again assume x > 0, as the 5-dimensional geometry 
is obviously singular at r = if x < 0. In the case x = 0, then y = and b = ±a, leading 



to the 5-dimensional "extreme black hole" metric with flat spatial sections [ft 13 



ds * = + dt 2 T — dtdx 5 - (l - - ) (dx 5 ) 2 - dr 2 - r 2 dtf . (3.11) 



While from (|3.5|) the Kretschmann scalar vanishes in this case, the study of geodesic 
motion shows that this geometry is actually singular. The effective potential 

v = t] + n 5 2 - n 4 2 + (114 ± n 5 ) 2 - (3. 12) 

r 

is unbounded from below for a < 0, in which case all the / = geodesies terminate at the 
singularity r = 0. For a > 0, only geodesies with / = n 4 ± IT5 = 0, rj < terminate at the 
singularity, which is thus harmless to physical test particles. 

The general case x > is quite similar to the case b) above, the 5-dimensional 
spacetime being geodesically complete, with an infinite number of changes of sign of A55. 



3.3 Case y > x 2 . 

Inspection of the 5-dimensional metric (|2.7|) , (|2.28|) with a given by (|2.25|) shows that for 
all values of x it is regular for r e] — 00, +00 [ , so that the 5-dimensional geometry 
is of the Lorentzian wormhole type. However, as pointed out in |18| (in the special case 
x = of the symmetrical wormhole), these wormholes are non-traversable, in the sense 
that physical (non-tachyonic) test particles cannot go from one asymptotically flat region 
(r — > +00) to the other (r — > —00). This would be possible only if the effective potential 
V(r) in ( |3.7|) was negative over the whole range of r. Here the potential ( |3.8| ) with 77 > 
is necessarily positive over part of this range, as the range ] —pTi/iu,, 0[ of pa is at least 2tt 
(p 2 > 4/i 2 ), and a physical test particle coming from r — > +00 is always reflected back to 
r — > +00, just as in the case x 2 /4 < y < x 2 . 
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Likewise, 5-dimensional light cones also gradually tumble over when r decreases from 
+00 to —00, leading to the existence of closed timelike curves. The 5-dimensional metric 
is asymptotically Minkowskian at both points at spatial infinity if p = 2np (y/x 2 = 
(n 2 — 1/4) /(n 2 — 1)) with n integer, in which case light cones tumble over by mr between 
the two points at infinity, and the solution is a metrical kink of winding number n |19|] 
(other axisymmetric kink solutions of Kaluza-Klein theory are discussed in ||20|| ). 



4 Small perturbations 

In this section we set up the equations for small time-dependent spherically symmetric 
perturbations of the electrostatic solutions. As shown in M, we may choose for the five- 
dimensional metric a parametrization similar to that of (2.2), ( ^29 ) with metric functions 
now depending on time, 

dsl = e 2 ^ dt 2 - e 2a ^ dp 2 - e 2 ^ dtf - e 2 ^ (dx 5 + 2V(j>, t) dt) 2 . (4.1) 

The Kaluza-Klein Gauss law (Eq. ( |2.4j )) may be integrated to 

F 41 = Q e- a - 2/3 -^ , (4.2) 

with Q the conserved electric charge. The remaining time-dependent spherically sym- 
metric equations ( |2.3| ) and ( |2.5| ) then reduce to the system 

e _ 2/3+ a- 7 -^ ( e 2/3-a+ 7 +^ ^y _ e 2(a- 7 ) ^ = _ 2 Q2 e 2( a -2/3-2^) ? (4 3) 

e -2/3+a- 7 -V^ e 2/3-a+y+i) p/y _ g2(a- 7 ) jj _ g2(a-/3) ^ (a ^\ 

e _ 2/3+a - 7 -^ e 2/3- a+7 +^ yy _ e 2(a-7) ( a + 2[3 + 4>) = 2Q 2 e^" 2 / 3 " 2 ^ , (4.5) 

-2(3" - if/' + (3' (2a' - 3/3') + ip'(a' - 2(3' - if,') + e 2 ^~^ = Q 2 e*(°>-W-*l>) , (4.6) 
-2$' {2(3' + if)')a + 2{-(3' + i)(3 + (Y - if)')^ = (4.7) 

with ' = d/dp, ' = d/dt. These equations are not all independent. The first three dy- 
namical equations correspond respectively to Eq. (|2.5|) and to the R\ and R\ components 



i, 2 aim 

of Eq. (|2.3| ), while the last two constraint equations correspond to (Rf — RD/2 and _R 4 i 
respectively. 

Now we linearize the metric fields in ( f4.1|) around the static background fields in Eqs. 
fl2T2"9] ), 

if;(p,t) = if;(p) + 6if;(p,t), etc, (4.8) 

where 5if> ,8a ,5(3 , £7 are small perturbations. The linearization of equations ( f4.3|) -( fO|) 
leads to the differential system for the perturbations, which has been simplified by using 
the static equations of motion, and choosing a "harmonic" background coordinate system 
in which the background fields if) , a , (3 , 7 are related by 

a- 2(3 - 7 - <0 = (4.9) 
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(which amounts to choosing the radial coordinate p = er, with h pp = H 2 ): 



Sip" - ip'5a' + 2ip'8f3' + ip'5i + ip'Sip' 

-e 2 W+^5i> - 2i\)"{5a - 25(3 - 2Sip) = , 
5(3" - (3' 5a' + 2(3' '5 (3' + (3'8i + (3'5i(j' 

-e 2 ^+^5(3-2(3"{5a-5(3)=0, 
5i' + i{-5a' + 25(3' + Si + Stf) - e^ 2f3+ ^(Sa 

+25(3 + Sip) + 2i("{5a - 25(3 - 25ip) = , 
25(3" + 5iP" - {2(3' + ip')Sa' + 2{(3' - i)5(3' + [ip' - i)Sip' 

-{2(3" + iP")5a + 2{(3" + i(")5(3 + 2iP"5iP = , 
25(3' + &tj>' = {2(3' + iP')Sa + 2{-(3' + i)5(3 + {i - ip')Sip . 



(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 



The last two constraint equations are not independent, as the time derivative of (|4.13|) 
may be seen (using 7" + ip" = 0, which follows from (|2.30|) ) to be identical with the space 
derivative of ( |4.14j ). 

The perturbations 5ip, 5a, etc. are defined only up to a change of coordinates pre- 
serving the form of the metric (|4.1| ), so that we still have the "gauge freedom" to choose 
coordinates for the perturbed spacetime by imposing a supplementary relation between 
these perturbations. The separation of the linearized equations ( 4. 10| )— ( 4.14 ) is simpler in 
the gauge where the perturbations 5\ip, 5ia, etc. are constrained by the gauge condition]] 



25t(3 + 5np = 0. 



The integration of the constraint equation ( 4.14 ) then leads to the relation 

f-(3' 



5\a 



ip' + 2(3' 



(4.15) 



(4.16) 



With the help of the unperturbed field equations and Eq. ( |4.15| ), Eqs. ( |4.10| ), ( ^4.11| ) 
combined together according to: (3'x ( |4.10D — ^ x (|4.11] ), lead to the wave equation for 

Slip 

S^" _ e 2 ^ + ^S^ -6^5^ = 0, (4.17) 



F 2 



with 



F = 2/ip' + 1/(3' . (4.18) 
^From a solution 5iip{a,t) to this master equation, the other perturbations in the gauge 



(|4.15| ) may be obtained by using Eqs. Q4.15| ), Q4.16| ), and the equation 



Sii 



ip' + 2(3' 



7 Sii/;' - <5ia' 



(4.19) 



obtained by adding Eqs. (|4.10|) and (|4.11| ) according to (|4.10|) + 2 x (|4.11|) , and using, 
first Eq. (|4.15|) to simplify the obtained equation, then Eq. ([4.16|) to cancel ip" , (3" . 



In the Einstein frame (J2.6|) this gauge condition reads simply 8\(3e = 0. 
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One can also gauge transform these perturbations 8ii(, Sio>, etc. to obtain the pertur- 
bations 5ip, 8a, etc. in a generic gauge, by carrying out the coordinate transformation 
(ti, pi) — ► (t, p). For instance the scalar field ip can be linearized around its static value 
in both coordinate systems, 



ij)(t,p) = ip(p) + 5ip(t, p)+-- ■ 

= i)(pi) + 5iip{ti,pi) H 

= ip(pi) + tfipi) Sxpfa, Pi) + Stp(t, p) + ■ 



(4.20) 



where we have linearized the coordinate transformation according to p ~ p\ + 8\p{ti, Pi)- 
To first order Eq. ( 4.20 ) leads to 



Sip = 8iip — ij)' Sip . 



(4.21) 



Eliminating 8\p between ([4.21 ) and a similar equation for the scalar perturbation 8(3, and 
using the gauge condition ( |4.15|) , we obtain the relation 



W-W=-(2^ + V')^- 



(4.22) 



^From this last relation it can be shown that our wave equation (|4.17|) is identical with 
the scalar wave equations (18) or (20) of || (see also 0) derived in the gauge 8<i(3 = 0, 
the relation between the Kaluza-Klein scalar perturbations, 



81^ 



2(3' 



-S2IP • 



2(3' + ip' 

leading to the identification with the notations of M, 

Snp = -^/2 = -fR/2. 

with/ = e" 2 ^ 2/?'/(2/3' + if/) . 



(4.23) 



(4.24) 



5 Stability 

We now address the question of stability of electrostatic solutions against radial per- 
turbations. Given stationary perturbations of the form 8iip(a,t) = Siijj(a) e int , etc., we 
search for non-trivial real solutions to the equation obtained from ( |4.17|) by assuming Q 
imaginary, Q = —ik (k > 0), 

W - (k 2 e 2 ^) + QF'/F 2 ) 8^ = 0, (5.1) 



where F is defined in ( |4.18 ). The existence of such solutions, satisfying some physical 



boundary conditions, means that any initially small perturbation will grow exponentially 
in time, thus the background solution is unstable. Conversely, the background solution is 
stable if all the eigenvalues Q are real. 
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The effective potential function U(a) — the coefficient of 8%ijj in the second term of Eq. 
(|5.1| ) — has double poles at the zeroes <jj of F, i.e. the roots of 2(3' + ip' = 0. Near such 
possible poles, the effective potential behaves as 

u(a) - 12 A /2 1 _ 2 

with /3| = /3'(<Ji), etc., from which follows the behaviour 

8^ ~ + C 2 (o- - (5.3) 

(J - CTj 

(Ci and C2 integration constants). So 5i^> generically has poles at a — Oi. However such 
poles turn out to be spurious, being induced by the gauge fixing 25i/3 + Slip = 0. Indeed, 



it follows from Eq. (|4.22 ) that perturbations 8ip and 8f3 which are regular in a generic 



gauge lead, in the gauge ( |4.15| ), to a 8iip with poles at <7j. Conversely these poles can be 



removed by transforming to another gauge. Similarly, from Eq. (f4.23| ) spurious poles also 
occur in the gauge 82P = at the zeroes of f3' [|10|, HQ 

These spurious poles being discarded, divergences of the perturbation 8.\ij) can only 
occur at the two ends of the range l a of a, and must be excluded by the choice of 
appropriate boundary conditions. We shall adopt here, as physically reasonable from the 
5-dimensional general-relativistic point of view, the boundary conditions previously used 
in scalar-tensor theories [ ID|| , which state that the relative perturbations of the background 



fields must be finite at the boundary. Since our functions 8%a, 8i(3, #17, 8\ip are relative 
perturbations of the background metric fields in (|2.2j), these boundary conditions read 



\8\a\ < 00 , \8\P\ < 00 , |<5i7| < 00 , \8iip\ < 00 , for a G l a (5.4) 

(the last of these boundary conditions was termed "strong" boundary condition in Ref. 
p0| ). The boundary conditions ( |5.4| ) lead after gauge transformation to similar boundary 
conditions in a generic gauge. 

Eq. ( |5.1D behaves at spatial infinity (r — > +00, a — > 0_) as 

S^'-^Sxip-O, (5.5) 
a 

for all values of the parameters (x, y, a), leading to the asymptotic behaviour 

5ii/j ~ a (cie fc/CT + c 2 e- k/a ) . (5.6) 
Our boundary conditions are satisfied by choosing C2 = 0. 



A number of possible behaviours of Eq. (5.1) can occur at the lower end a m i n of the 



range I a , depending on the values of the parameters (x, y, a). In most cases 0, the 
corresponding general asymptotic solution will be a linear combination of a bounded and 



7 In Refs. [|| and M, solutions for which such poles occurred were argued to be stable, on the basis 
that generic perturbations S21P were unbounded; however this argument is not gauge invariant. 
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an unbounded solution with coefficients C3 and C4. Then the perturbation will remain 
bounded near cr min provided c 4 is fixed equal to zero. However, owing to the scale invari- 
ance of Eq. (|5.1|), the two conditions c 2 = c 4 = can be satisfied simultaneously only for 
special values of the parameter k 2 , i.e. we have an eigenvalue problem. In the absence 
of the knowledge of the exact solution in the whole interval, it is impossible in general to 
solve analytically this eigenvalue problem, and thus to establish with certainty the sta- 
ble or unstable character of the background solution. However there are two parameter 
ranges, corresponding to neutral solutions, and to a preferred set of charged solutions, for 
which stability can be proved analytically by straightforward arguments. We will examine 
these cases in the next two subsections, and discuss in a third subsection two special cases 
for which unstability can be proved analytically. 



5.1 Stability of neutral solutions 



In the parameter space (x, y, a), the domain of neutral solutions is obtained from (|2.21 
by setting 6 = 0. Since a is real, this leads to the portion of the surface 

a 2 — xa + y = , (5.7) 

for which y < x 2 /A (q 2 > 0, with q = \x/2 — a\). As recalled in the Introduction, neutral 
Kaluza-Klein theory corresponds to Brans-Dicke theory with u = 0. Therefore the 
proof, given in |TI| that all static spherically symmetric solutions to Brans-Dicke theory 
are stable under radial pertutbations carries over to the present case. For completeness, 
we sketch here this proof in the Kaluza-Klein setting. 

In the neutral case, the static scalar equation (|4.3| ) reduces to ip" = in the coordinate 
system ( f4.9p , so that the linearized equation ( |4.10| ) immediately decouples from the other 



equations in the harmonic gauge 

5 a - 25 p - 5 7 - 5 ip = . (5.8) 
The resulting wave equation for the perturbation S Q ip is 

W -k 2 e 2 ^ + ^5 ^ = 0. (5.9) 
The ratio 5qiJj" = k 2 e 2 ^ + ^ being regular and positive over I a =] — 00, 0[, it is 



impossible to keep finite at both ends of I a . Our boundary conditions (|5.4|) cannot 
be satisfied, so that this case is stable. 

Tomimatsu investigated the stability of 2-static solutions with vanishing electric 
field ([5.71). In the parameter space (x, y, a), he considered, as shown in M, only the surface 



portions: a = x/2 — q, y < 0, x > 0, and a = x/2 + q, y > , x < 0, with y < x 2 /4. 
He concluded that all the solutions under consideration were unstable, except the case 
y = a = 0, 1 ^ of the Schwarzschild solution which was found to be stable. The 
discrepancy with our conclusions can be explained, as discussed in [jj, by the fact that 
Tomimatsu treated the Schrodinger eigenvalue problem for an auxiliary function whose 
direct physical meaning is not transparent, and used in effect boundary conditions less 
stringent than our conditions ( |5.4| ). 
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5.2 A set of stable charged solutions 



Let us return to the wave equation ( |5.1|) in the gauge ( |4.15| ). The ratio is 
positive definite over l a (V/c 2 > 0) provided the two conditions 



,2V 



> 



F' > 



(5.10) 
(5.11) 



are satisfied simultaneously. Then, if the range l a is infinite (which is the case for y < x 2 ) 
and if further F does not have a zero in this range, 



F^O 



(5.12) 



it is impossible to satisfy our boundary conditions ( |5.4j ), and the corresponding static 
solution is stable (a zero of F would lead, as discussed above, to a pole of the perturbation 
5ii/j(a), which would enable it to remain finite at both ends of I a with S^'/Siip > 
everywhere; as we show in the Appendix the condition ( |5 . 1 2 ) is always satisfied when the 
conditions ( |5.10| ) and ( |5.11| ) are fulfilled). 

Eq. (ET28D shows that e 2 ^ = -A 55 chan ges sign periodically for y > x 2 /A, so that 
condition ( |5.10| ) cannot be satisfied, while it is satisfied for y < x 2 /4 if x/2 — a > 0, 
implying (on account of (|2.21|) ) 

x/2-a>\b\. (5.13) 
To investigate condition ( |5.11|) , we compute (for y < x 2 /4) 



P' 



X Q 

- + - coth(ga - rj) 



x 



v coth vo , 

2 



(5.14) 
(5.15) 



where i] is defined by q = \b\ sinhrj, x/2 — a 

,2 



|&| cosh?7 (rj > 0). Eq. ( [4.18| ) then leads to 



F' 



q 



ijj' 2 sinh (qa — rf) (3' 2 sinh (i/<t) 



(5.16) 



The determination of the range for which this function remains non-negative is carried out 
analytically in the Appendix. The conclusion is that our stability conditions ( 5.10| )-( |5.12| ) 
are satified if either 

a) x = a — b 2 /a (implying y = 0) with —2\b\ < a < 0. This class of stable solutions, 
includes black holes (x > 0, a < 0) and extreme black holes (x = 0, a < 0); 

b) 2(a + < x < a - b 2 /a (implying y > 0) with —2\b\ < a < —\b\. 



5.3 Charged solutions: two unstable cases, and a stable subcase 
5.3.1 Case y > 0, x = a = 

In this massless case (M = S = 0), the matrix A in (|2.28|) reduces to 



A 



cos bo 
sin bo 



sin bo 

- cos bo 



(5.17) 



15 



leading to a symmetrical Lorentzian wormhole spacetime, with 



,2a 



1, e 2f3 = r 2 + fi 2 , e^ = e- 2 ^ = (r 2 + /i 2 )/(r 2 - yU 2 ) 



(5.18) 



(r = — fx cot(/xcr), \x = b/2 = Q). A deeper investigation of this solution is given in JIB 



Due to the symmetry of the solution under the change in the radial coordinate r — > —r 
and — consequently — of the wave equation ( |5.1| ) (rewritten in terms of the coordinate r) 



[(r 2 + iJ 2 )S 1 i) >r ] )r - 



k (r 2 -/i 2 ) + 



6/i 2 



5^ = 



(5.19) 



the two end-points at spatial infinity (r — > ±oo) are merged, leading to only one divergence 
at r = +oo, which can be cancelled by choosing C2 = in Eq. (|5.6|) . At the lower end of 
the range [0, +oo[ of r, the perturbation 5iip is unbounded (cf. Eq. (29) of || and Eq. 

(EH)) 



Slip 



c 3 



i + ; k 2 + 



r 2 + 



+ c 4 r 6 



1 - 



14 



k - 



34 



r 2 + ■ 



(5.20) 



Actually the pole in (|5.20| ) is gauge dependent and can be removed by transforming to 
the gauge <5 2 /3 = 0, according to Eq. (|4.23|) , which leads to 



5 2 ip 



r 2 — ji 2 



Sup 



C3 
2 



i + ; k 2 + 



r 2 + ■ 



C4 r 



1 - 



14 



k - 



76 



r 2 + 



(5.21) 



This is well bounded for all values of the integration constants. However the constraint 
equation ( |4.14| ) now leads to 



5 2 a 



S 2 ^' + (V - Y)^ 



c 3 /j 
6r 2 



2^' + V>' 



which is singular, unless c 3 (/c) = (leading to an eigenvalue problem), or 

k = 2/(JL. 



(5.22) 



(5.23) 



We verify that — for this special eigenvalue — all the perturbations 5 2 ip, etc, are bounded 
Vr G [0, +00 [. The massless charged symmetric wormhole solution is then unstable for 
it admits a mode of perturbation growing in time as e 2 *^. 
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5.3.2 Case y — 0, a > x > 
The functions which appear in Eq. 

2lp 



are given by 



x 



1 

X 



X 



Ax 



a 

■ + -e 

x 

3a 



ax 



e~ xr7 + - 
x 



(5.24) 
(5.25) 
(5.26) 



With the help of Eq. (|A.2|) , we see that the third term in Eq. (pTT ) vanishes identically 
if 4x — 3a = 0, or behaves at the lower end of I a (a — > — oo) as: 



QF' 



6ax 2 



Ax 



F 2 3a 

if Ax — 3a ^ 0. Since the second term in Eq. (|5.1| ) behaves as a constant (a —>■ — oo): 

e 4/3+2 ^~-x 3 (a-x), 
we can neglect the third term, and write Eq. ( |5.1D as 

5^" + k 2 x 3 (a - x) 5^ ~ . (5.27) 

We obtain by integration 

5i^> ~ C3 cos ma + C4 sin ma (5.28) 

(with m = kx\jx[a — x)), which is bounded for all values of the integration constants. 
We must check that the other perturbations 8\a and 8yy are also bounded. Using the 
asymptotic behaviour 



i)' + 2(3' 



x — 2a 
3a — Ax 



1 

- e 
2 



if 3a - Ax ^ 



if 3a — Ax = . 



(5.29) 



we find from Eqs. (|4.16 ) and ( |4.19|) that for a ^ Ax/ 3 all the perturbations are finite at 
a — > —00. It follows that the boundary conditions ( |5.4D are satisfied for all values of k if 
c 2 = in ( |5.6| ), and the corresponding background solutions are unstable. 

However if a = Ax/3 = 2\b\ (E = -2y/5M, \Q\ = 2M), 5 x a and 5^ diverge, from 
( 4.16 ), ( 4.19 ) and ( |5.29| ), as e~ xa 5\ip. These divergences may not be removed by trans- 
forming to another gauge. For instance, in the gauge 5 2 f3 = 0, S^ip ~ e xa Slip from ( [4.23|) , 
leading from the first equation ( |5.22| ) to a 5 2 a diverging as e~ xa times a bounded function. 
A similar divergence results in a generic linear gauge {8f3 + \8ip — 0) from the analysis of 
Eqs. ( |4.14j ) and ( |4.22j ). So at least one of the perturbations is never bounded^, leading to 
the conclusion that this special subcase is stable. 

8 One can show that 5a is bounded in the gauge 5a — Sf3 = 0; however Sj is not bounded in this gauge. 
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6 Conclusion 



We have discussed the geometry of the 3-parameter family of static spherically symmetric 
solutions of 5-dimensional Kaluza-Klein theory. We have found that, besides the 2- 
parameter black hole class (y = 0, x > 0, a < 0) and the exceptional regular solution 
(y = 0, a = x > 0), this family contains a 3-parameter class of geodesically complete 
solutions (y > 3x 2 /4, x > 0), which is larger than previously thought. These regular 
solutions are necessarily charged. 

We have also investigated the stability of these solutions under radial perturbations. 
We have shown that all neutral solutions (b = 0) are stable. Among charged solutions, 
we have found two stability classes (y = 0, a < 0, x > 3a/ 4) and (0 < y < x 2 /4, —2\b\ < 
a < —\b\); the first stablity class includes all the black hole and extreme black hole 
solutions. We have also been able to prove analytically the unstability of two lower- 
dimensional (in the 3-dimensional parameter space) classes of solutions. The remaining 
cases (which include the geodesically complete solutions) lead to eigenvalue problems. 
Pending a numerical investigation of these, we can only conjecture that the corresponding 
static solutions are unstable. 

This work should be generalized in two directions. Electrostatic spherically symmetric 
solutions are only a subcase of the stationary spherically symmetric solutions which may 
be obtained from Eq. ( |2.12| ), and which are discussed in the extreme case in [|13|, II |. The 



stability of these solutions could be investigated along the lines followed here. Also, the 
proof of stability under monopole perturbations is but a first step. The effect of higher 
multipole perturbations should also be ascertained before one can conclude to stability 
under all small perturbations, but this is a far more difficult task. 
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Appendix 



In this Appendix we determine (for y < x 2 /4 and x/2 — a > \b\) the parameter domains 
for which the function 

Q 2 V 2 

^ ^ -ip' 2 sinh 2 (ga — rj) (3' 2 sinh 2 (z/a) ^ ^ 

remains non-negative in the range ] — oo, 0]. There are 3 cases according to the value of 

y- 
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1) y < (g > v) . For a — > — oo, F' ~ — u 2 / (3' 2 sinh 2 z/cr, so that condition (|5.11 ) is 
not satisfied. 

2) y = (q = v = \x\/2), i.e. 5 2 = a(a — x). Then 

F' = 3a ~ AX e- x °. (A.2) 



Remembering that a < x/2, we find that condition Q5.ll ) is satisfied for a < 0, x > 3a/4 



implying 6 2 > a 2 /4. Condition ( |5.12| ) is then also satisfied because 

F{a) < F(0) = 2/<0'(O) = 4/o < . (A.3) 
The corresponding domain of stability is 

x = a-6 2 /a, -2|6|<a<0. (A.4) 
3) y > (q < v). Then, Eq. ( A.l ) can be written as 



4 sinh 2 5 sinh 2 k ^. . 

sinh 2 (gcr — rj + e5) sinh 2 (z/cr + en) 

where 5 > and k > are defined by 

x 

2 =£ v / ^ cosh5 ' g = v/ysinh(5 

— = cosh k , i/ = sinh k 
2 2' 2 

(e = sign(x)). If x > 0, ua + K has a zero for a = — k/V, with F' — >■ — oo, so that condition 
(|5.11| ) is not satified. 

If x < 0, ^'(0) = (46 2 — a 2 ) /a 2 is non-negative if 

6 2 >a 2 /4. (A.6) 

F'(a) can vanish only at the zeroes of 

h = f-g, (A.7) 

where the functions 

/ = — sinh(z/cr — k) , <7 = sinh(go" — r] — 5) . (A. 8) 

g 

are positive for a negative. However h(a) stays positive in l a because the inequalities 

(u - q)a < < k - r] - 5 (A.9) 

(see below) imply 

h' = — z/[cosh(z/er — k) — cosh(ga — rj — 5)] < , (A. 10) 
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so that the minimum value of h is h(0) = (i//^/y)(2 + a/|6|) > 0. Therefore F' stays 
positive in I a , so that (condition (|5.12| ) being again satisfied by virtue of (|A.3| )) the 
corresponding static solutions, with parameters 



2(a+ |6|) < x < a-b 2 /a, -2\b\ < a < -\b\ 



(A.ll) 



(the first two inequalities result from ( |5.13| ) and y > 0, the last two from (A..6) and x < 
together with (|5.13|) ) are stable. 

It remains to prove the second inequality ( |A.9| ). From the definitions of rj, k and 5 we 
obtain 



\b\y sinh(/t — 77 — 5) 



= (u — q)ax — 2vy 

= 2u(a 2 — b 2 ) — {y + q)ax 

> 2(|o|-|6|)(i/|6|-g|o|) ) 



(A.12) 
(A.13) 



using the first inequality ([A.ll ). The first factor is non-negative from the last inequality 
(lATHD , while 



2r,2 2 2 

v b — q a 



(x 2 /A){b 2 - a 2 ) + y{a 2 - 6 2 /4) > (36 2 /4)y > 



(A.14) 



for < y < x 2 /4. 
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